figure(53)
FigureTitle = 'Figure11-CrustalThickness-Trace';
set(gcf,'name',regexprep(FigureTitle,'\_*','-'))


            hold on
            subaxis(3,2,1)
            hold on
            axis([0 11 .6 .8])
            ax = gca;
            set(gca, 'Xtick',[0:1:11])
            ax.XAxis.MinorTick = 'on';
            increment = .5;
            ax.XAxis.MinorTickValues = ax.XLim(1):increment:ax.XLim(2);
            set(gca,'Ytick',[0:.05:1])
            ax.YAxis.MinorTick = 'on';
            increment = .01;
            ax.YAxis.MinorTickValues = ax.YLim(1):increment:ax.YLim(2);
            
            errorbar(6.5, 0.775, .456, .456, 1, 1,'LineWidth',5,'Color','k')
            
            
%             ylabel(sprintf('%s_{%s}',RatioLabels{1}, normalizationLabel));
%             yline(0.775,'k-','N-MORB','LineWidth',2);
%             yline(0.63,'k-','DMM','LineWidth',2);
%                         xline(6.5,'k-','LineWidth',2);
%             xline(5.5,'k-','LineWidth',1);
%             xline(7.5,'k-','LineWidth',1);
%             yline(1,'k:','PUM','LineWidth',2);
           
            subaxis(3,2,2)
            hold on
            axis([0 11 .6 1.4])
            set(gca, 'Xtick',[0:1:11])
            set(gca,'Ytick',[0.5:.1:2])

            ax = gca;
            ax.XAxis.MinorTick = 'on';
            increment = .5;
            ax.XAxis.MinorTickValues = ax.XLim(1):increment:ax.XLim(2);
     
            ax.YAxis.MinorTick = 'on';
            increment = .1;
            ax.YAxis.MinorTickValues = ax.YLim(1):increment:ax.YLim(2);
            
            errorbar(6.5, 1.153, .044, .044, 1, 1,'LineWidth',5,'Color','k')
            
%             ylabel(sprintf('%s_{%s}',RatioLabels{11}, normalizationLabel));
%             yline(1.152,'k-','N-MORB','LineWidth',2);
%             yline(0.71,'k-','DMM','LineWidth',2);
%             xline(6.5,'k-','LineWidth',2);
%             xline(5.5,'k-','LineWidth',1);
%             xline(7.5,'k-','LineWidth',1);
%              yline(1,'k:','PUM','LineWidth',2);

             yline(0.71,'k-','WHDMM','LineWidth',2,'Color','k','FontName','Times New Roman');
             yline(0.668,'k-','LineWidth',1,'Color','k');
             yline(0.744,'k-','LineWidth',1,'Color','k');
            

            yline(0.937,'k:','MK91','LineWidth',1,'Color','k','FontName','Times New Roman');
            yline(0.798,'k:','HS96','LineWidth',1,'Color','k','FontName','Times New Roman');
            yline(0.851,'k-','This Study','LineWidth',3,'Color','k','FontName','Times New Roman');

            
            
            

            subaxis(3,2,3)
            hold on
            axis([0 11 .7 1.7])
            ax = gca;
              set(gca, 'Xtick',[0:1:11])
            set(gca,'Ytick',[0.5:.1:2])

            ax = gca;
            ax.XAxis.MinorTick = 'on';
            increment = .5;
            ax.XAxis.MinorTickValues = ax.XLim(1):increment:ax.XLim(2);
     
            ax.YAxis.MinorTick = 'on';
            increment = .1;
            ax.YAxis.MinorTickValues = ax.YLim(1):increment:ax.YLim(2);
            
            errorbar(6.5,0.813, .011, .011, 1, 1,'LineWidth',5,'Color','k')
                        
%             yline(1,'k:','PUM','LineWidth',2);
%             yline(0.82,'k-','N-MORB','LineWidth',2);
%             yline(1.5,'k-','DMM','LineWidth',2);
%             
%             xline(6.5,'k-','LineWidth',2);
%             xline(5.5,'k-','LineWidth',1);
%             xline(7.5,'k-','LineWidth',1);
%             yline(1,'k:','PUM','LineWidth',2);


            yline(1.53,'k-','WHDMM','LineWidth',2,'Color','k','FontName','Times New Roman');
            yline(1.46,'k-','LineWidth',1,'Color','k','FontName','Times New Roman');
            yline(1.65,'k-','LineWidth',1,'Color','k','FontName','Times New Roman');
            
            
            yline(1.03,'k:','MK91','LineWidth',1,'Color','k','FontName','Times New Roman');
            yline(1.28,'k-','HS96,This Study','LineWidth',3,'Color','k','FontName','Times New Roman');
            yline(1.18,'k:','CBT01','LineWidth',1,'Color','k','FontName','Times New Roman');

            
            

            subaxis(3,2,4)
            hold on
            %axis([0 .15 1 4])
            axis([0 11 .45 .9])
            ax = gca;
              set(gca, 'Xtick',[0:1:11])
            set(gca,'Ytick',[0.45:.05:.9])

            ax = gca;
            ax.XAxis.MinorTick = 'on';
            increment = .5;
            ax.XAxis.MinorTickValues = ax.XLim(1):increment:ax.XLim(2);
     
            ax.YAxis.MinorTick = 'on';
            increment = .01;
            ax.YAxis.MinorTickValues = ax.YLim(1):increment:ax.YLim(2);
            
            
            errorbar(6.5,0.753, .124, .124, 1, 1,'LineWidth',5,'Color','k')
                
%             yline(0.59,'k-','N-MORB','LineWidth',2);
%             yline(0.5,'k-','DMM','LineWidth',2);
%             
%                         xline(6.5,'k-','LineWidth',2);
%             xline(5.5,'k-','LineWidth',1);
%             xline(7.5,'k-','LineWidth',1);
%             yline(1,'k:','PUM','LineWidth',2);
            

            yline(0.503,'k-','WHDMM','LineWidth',2,'Color','k','FontName','Times New Roman');
            yline(0.45,'k-','LineWidth',1,'Color','k','FontName','Times New Roman');
            yline(0.54,'k-','LineWidth',1,'Color','k','FontName','Times New Roman');

            yline(0.43,'k:','MK91','LineWidth',1,'Color','k','FontName','Times New Roman');
            yline(0.62,'k-','This Study','LineWidth',3,'Color','k','FontName','Times New Roman');
            %yline(1.18,'k:','CBT01','LineWidth',1,'Color','k','FontName','Times New Roman');
            
            
            
            
            
            subaxis(3,2,5)
            hold on
            axis([.1 2 .5 1.5])
            set(gca, 'Xtick',[0.1:.2:2])
            set(gca,'Ytick',[0.5:.1:2])
            ax = gca;
            ax.XAxis.MinorTick = 'on';
            increment = .1;
            ax.XAxis.MinorTickValues = ax.XLim(1):increment:ax.XLim(2);
     
            ax.YAxis.MinorTick = 'on';
            increment = .1;
            ax.YAxis.MinorTickValues = ax.YLim(1):increment:ax.YLim(2);
            
         
            errorbar(0.813,1.004, .018, .018, .011, .011,'LineWidth',5,'Color','k')
%             xline(1,'k:','PUM','LineWidth',2);
%             xline(0.82,'k-','N-MORB','LineWidth',2);
%             xline(1.548,'k-','DMM','LineWidth',2);
%             
%             yline(1,'k:','PUM','LineWidth',2);
%             yline(1,'k-','N-MORB','LineWidth',2);
%             yline(1.3,'k-','DMM','LineWidth',2);
            

            xline(1.53,'k-','WHDMM','LineWidth',2,'Color','k','FontName','Times New Roman');
            xline(1.46,'k-','LineWidth',1,'Color','k','FontName','Times New Roman');
            xline(1.65,'k-','LineWidth',1,'Color','k','FontName','Times New Roman');
            xline(1.03,'k:','MK91','LineWidth',1,'Color','k','FontName','Times New Roman');
            xline(1.28,'k-','HS96,This Study','LineWidth',3,'Color','k','FontName','Times New Roman');
            xline(1.18,'k:','CBT01','LineWidth',1,'Color','k','FontName','Times New Roman');

 
            yline(1.27,'k-','WHDMM','LineWidth',2,'Color','k','FontName','Times New Roman');
            yline(1.22,'k-','LineWidth',1,'Color','k','FontName','Times New Roman');
            yline(1.32,'k-','LineWidth',1,'Color','k','FontName','Times New Roman');
            yline(1.21,'k:','HS96','LineWidth',1,'Color','k','FontName','Times New Roman');
            yline(1.13,'k-','MK91, CBT01, This Study','LineWidth',3,'Color','k','FontName','Times New Roman');
            
            
            
            
            subaxis(3,2,6)
            hold on
            axis([0 1 0 .6])
%             xlabel(sprintf('%s %s_{%s}','\delta',xxvalue{:}, normalizationLabel));
%             ylabel(sprintf('%s %s_{%s}','\delta',yyvalue{:}, normalizationLabel));
            ax = gca;
            set(gca,'Xtick',[0:.2:1])
            ax.XAxis.MinorTick = 'on';
            increment = .1;
            ax.XAxis.MinorTickValues = ax.XLim(1):increment:ax.XLim(2);
            
            set(gca,'Ytick',[0:.1:1])
            ax.YAxis.MinorTick = 'on';
            increment = .05;
            ax.YAxis.MinorTickValues = ax.YLim(1):increment:ax.YLim(2);
            
            %axis([0.05 0.5 0 .5])
%             ax = gca;
%             %increment = findBestIncrement( ax.YLim(2)- ax.YLim(1),ax.YAxis.TickValues(2) -  ax.YAxis.TickValues(1));
%             set(gca,'Ytick',[.2 .5 1 1.3 2 3])
%             ax.YAxis.MinorTickValues = [.2 1 3];
%             ax = gca;
%             set(gca,'Xtick',[0:2:20])
%             ax.XAxis.MinorTick = 'on';
%             increment = 1;
%             ax.XAxis.MinorTickValues = ax.XLim(1):increment:ax.XLim(2);
%             yline(1,'k-','Secular Equil','LabelHorizontalAlignment','right','LineWidth',2);    
%                         yline(1.2,'k:');
%             yline(1.25,'k:');
%             yline(1.13,'k:');
            

            for kk=1:6
                subaxis(3,2,kk)
                hold on
                ax=gca;
                set(gca,'ticklength',1.5*[0.0300    0.0500])
                set(gca,'fontsize', 15,'LineWidth',0.7,'FontName','Times New Roman')
                %set(gca,'yscale','log')
                 ax.YAxis.MinorTick = 'on';
%                ax.YAxis.MinorTickValues = [0.001:0.001:.01 .02:.01:.1 .2:.1:1 2:1:10 20:10:100];
%                 curtick = get(gca, 'XTick');
%                 set(gca, 'XTickLabel', cellstr(num2str(curtick(:))));
%                 curtick = get(gca, 'YTick');
            %    set(gca, 'YTickLabel', cellstr(num2str(curtick(:))));
                set(gca,'XColor', 'k')
                set(gca,'YColor', 'k')
            end

                 